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Abstract 

In this work we analyze the emission of gravitational waves from a gravita- 
tional system described by a Newtonian term plus a Henon-Heiles term. The 
main concern is to analyze how the inclusion of the Newtonian term changes 
the emission of gravitational waves, considering its emission in the chaotic 
and regular regime. 
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In a previous paper we analyzed the emission of gravitational waves from a Henon- 
Heiles System, showing the qualitative differences between gravitational waves emission 
from chaotic and regular system ]]]. Although interesting, the Henon-Heiles system is not 
realistic and it was originally intended to study galactic dynamics, such that the emission of 
gravitational waves is expected to be negligible. Thus in this work we include a Newtonian 
term oc r _1 to the Henon-Heiles potential. This class of potential was analysed by Vieira 
and Letelier ||, and correspond to a black- hole with external halo. This type of system is 
more realistic than simple Henon-Heiles system and its gravitational waves emission may be 
important due to its stellar system size. 

In this work, the dynamics are determined by the gravitational potential V = 
—GMm/r + mu 2 [(x 2 + y 2 )/2 + [x 2 y — y 3 /3)/a\. Using the usual approach, we work in 
a dimensionless system, such that the gravitational potential is reduced to 



$ = -— 

r 



(i) 



2 ' * 3 

with $ = V/(mLJ 2 a 2 ), a = GM/(u 2 a 3 ),r = r/a ,rj = x/a and £ = y/a. A main difference 
between pure Henon-Heiles 0] and above potential, is that bounded trajectories in the 
first case are limited to the range < E < 1/6, whereas in our case we have bounded 
trajectories even with E < 0. Other difference is that the upper limit in the energy ( 
E m ax = 1/6 ~ 0.167) may be lower than a pure Henon-Heiles case, its value depending of 
the a parameter. For example with a = 0.1, we have E max ~ 0.07 and with a = 0.01 we 
have E max ~ 0.156. With above potential, we calculated the trajectories using numerical 
simulations, changing energy E, the parameter a and initial conditions. For our simulations 
we used a = 0.01,0.001 and 0.0001. The main reason to these choice is to maintain some 
resemblance with Henon-Heiles case. (Later in the text, we discuss the physical meaning 
and limitations of this choice.) Note that with these a we have always E max ~ 1/6. We 
used in our simulations three distinct values for the energies, i.e., E = 0.1,0.01 and 0.001. 
In this range of energy we obtain a good mixture of chaotic and regular motions for our 
choice of as, permitting a better comparison between emission of gravitational waves from 



both cases. 

The power emitted as a gravitational waves was determined using the standard 
quadrupole formula 



dE , -2 



where 

Gm 9 n 

45c 5 v ' 

and quadrupole terms are given by q m = 2rf — £ 2 , = 2£ 2 — rf and q n % = 3r?£. (In what 
follow we use T instead of ui, with T = 2ir/u.) 

In the figures |l|J^ and |3| we show some results of our simulations. We plot trajectories 
in the X-Y plane, Poincare section and the spectrum of the emitted gravitational waves for 
a = 0.01 and three different energies E = 0.10,0.01 and 0.001. In each case we plot one 
chaotic and one regular motion. The maximum power emitted as a gravitational waves in 
each case are given in the table |. Note that in all cases, except when a = 0.01 and E = 0.1, 
the power emitted in chaotic case is greater than regular case ( This is an exception, as can 
be seen in the table [TV] where we show a more detailed statistics of emitted power.). 

The irregular characteristic of the gravitational waves spectrum is show in all chaotic 
case, being more strong with E — 0.1. In particular note the strong differences of the emitted 
spectrum of the chaotic and regular case in the figure |l|. 

In the tables ||,|T| and |V| we show with more details the power emitted as a gravitational 



waves for several a and E. As expected, the emission in the chaotic case is greater than 
regular case [|TJ. However, for lower E the differences between regular and chaotic case 
becomes negligible. This happens because for high E, the trajectories reach the boundaries 
of the Henon-Heiles term. In this case the trajectories of the particles suffer a strong changes, 
resulting in a higher emission of gravitational waves. On the other hand for lower E, the 
trajectories are confined far away from Henon-Heiles boundaries, such that it not suffer 
strong changes as the previous case. However, in all cases the spectrum of the emitted 



gravitational waves are much more irregular in the chaotic case when compared with the 
regular case. 

As soon the potential which we are using is more realistic than a simple Henon-Heiles case, 
it will be interesting to verify in which conditions the emission of gravitational waves will 
be important. To do it, it is important to fix the values of a and A. Thus in order to obtain 
a better understanding of a and A, we assume two typical size for the system: (1) Solar 
System and (2) binary pulsar (typically PSR 1913+16) with parameters (a±, Ti), (a 2 , T 2 ) 
respectively. Then we set T = T Q Tj , a = a cii , i = 1, 2 with a\ = 1.5 x 10 n m, T\ = 3.1 x 10 7 s, 
a 2 = 2.0 x 10 9 m and T 2 = 2.8 x 10 4 m and M = 3 X M Q , m = 3 2 M Q with M = 2.0 x 10 30 # 
(Solar mass). Thus 

T 2 a 2 
Solar System size : a x = 0.96/3iHr ! Ai = 2.29 x 10~ 22 3 2 -^ (4) 

T 2 a 2 
Binary Sistem size : a 2 = 0.33/3iHr ; A 2 = 5.5 x 10" 17 /5 2 -^ (5) 

In the tables [V] and [VT] we show some numerical values for a a , T Q and A using a = 
0.0001, 0.001 and 0.01 and 3 X = 3 2 = 1.4. Also, we show AE/Ei ~ (< P rnax > T a )/Ei which 
is the energy lost as a gravitational waves during a typical crossing time T and we use two 
values for Ei, i.e., Ei = E max = 0.10 and Ei = E m i n = 0.001. This last column show how 
much the emission are important to the dynamics of the system. Only when AE/E 1 the 
emission of gravitational waves will be negligible to the dynamics of the system. (Although 
AE/E Q > 1 is meaningless, we show these values only for the sake of completeness.) From 
the tables and we see that when the system size is of the order of a binary pulsar the 



emission of gravitational waves will be important. On the other hand when the size is of 
the order of the Solar System the energy lost as a gravitational waves are negligible only in 
the cases (la) and (lb). In all other cases the energy lost as a gravitational waves are of the 
same order of magnitude than its initial energy. 

Note that AE/E is lower when a and E are lower. As said before, in these cases the 
size of typical trajectory is small compared with the size of external halo. For example when 
a = 0.0001 and energies E D = 0.001 and E = 0.10 we obtained with numerical simulations 



that |f | 

max ^ 

0.07 , and |f | 

max ~ 0.6 r6sp6ctively. Thes6 result show that the size of the 
trajectories with small energy are smaller than trajectories with high energy (as expected). 
This mean that in the first case (small energy) the trajectories do not reach the boundaries 
defined by the Henon-Heiles term, whereas in the high energy case these boundaries can 
be reached. Thus in the first case (low energy), the external halo is felt only as a small 
perturbations (note that for small (rj,£) the Henon-Heiles term is negligible in respect to 
the Newtonian term). 

Therefore a potential like as given by equation [l], is useful only when external halo is felt 
as a small perturbations (as soon we consider the effect of the gravitational waves to the 
dynamics of the system). 

However, even with inclusion of the Newtonian term, the emission of gravitational waves 
are greater in the chaotic case than regular case - except for low energy case - , like our 
previous results 0. Also the spectrum of the emitted wave is highly irregular, thus we can 
speculate that with a careful analyze of the gravitational waves spectrum we can detect a 
gravitational chaotic system. But it is important to note that our results show that if a 
system is described by a Newtonian potential with a Henon-Heiles term its typical size need 
to be greater than a Solar System size. If this is not the case, the emission of gravitational 
waves will be important and it will collapse quickly. This show that although more realistic, 
a single inclusion of the Newtonian term is not sufficient as a realistic core-halo system model 
(see for example || for a more realistic approach). 

This work was partially supported by Fundagao de Amparo a Pesquisa do Rio Grande 
do Sul (FAPERGS), brazilian financial agency. 
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FIGURES 

FIG. 1. An examples of chaotic and regular motions with E = 1CT 1 and a = 1CT 2 . The left 
box show the trajectories in the X-Y plane, the central box the Poincare section and the right box 
the some portion of the gravitaional waves spectrum. The chaotic case is show in the top figures 
set, and regular case at bottom set. 

FIG. 2. Same as previous figure, but now using E = 10 -2 . 

FIG. 3. Same as previous figure, but now using E = 10~ 3 . 
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TABLES 



E Prnax / A 

0.100 5.21 x 10 17 chaotic 

0.100 5.34 x 10 17 regular 

0.010 1.49 x 10 16 chaotic 

0.010 6.73 x 10 15 regular 

0.001 1.04 x 10 16 chaotic 

0.001 5.23 x 10 15 regular 
TABLE I. Maximum emiited power as a gravitational waves for the cases shown in the figures 
21 and H 



a E P max /A Note 

0.0001 0.100 (8.15 ± 3.78) x 10 12 chaotic+regular 

0.0001 0.100 (5.89 ± 2.49) x 10 12 regular 

0.0001 0.100 (1.00 ± 0.36) x 10 13 chaotic 

0.0001 0.010 (6.26 ± 2.74) x 10 10 chaotic+regular 

0.0001 0.010 (6.43 ± 2.81) x 10 10 regular 

0.0001 0.010 (5.65 ± 2.43) x 10 10 chaotic 

0.0001 0.001 (2.72 ± 1.22) x 10 9 regular 

TABLE II. Maximum power emitted as a gravitationala waves with a = 0.0001 for three 
energies E = 0.1,0.01 and 0.001 



S 



a 


E 


Pmax 1 A 


Note 


0.0010 


0.100 


(1.53 ±0.73) x 10 17 


chaotic±regular 


0.0010 


0.100 


(9.61 ± 5.44) x 10 16 


regular 


0.0010 


0.100 


(1.78 ±0.66) x 10 17 


chaotic 


0.0010 


0.010 


(2.09 ±0.96) x 10 15 


chaotic±regular 


0.0010 


0.010 


(1.25 ±0.45) x 10 15 


regular 


0.0010 


0.010 


(2.31 ±0.94) x 10 15 


chaotic 


0.0010 


0.001 


(4.97 ± 2.41) x 10 14 


chaotic±regular 


0.0010 


0.001 


(5.00 ± 2.60) x 10 14 


regular 


0.0010 


0.001 


(4.87 ± 1.70) x 10 14 


chaotic 


TABLE III. Maximum power emitted 


as a gravitationala waves with a 


= 0.001 for three 


energies E = 


0.1,0.01 and 0.001 






a 


E 


Pmax 1 A 


Note 


0.0100 


0.100 


(2.55 ±0.79) x 10 17 


chaotic±regular 


0.0100 


0.100 


(1.41 ±0.59) x 10 17 


regular 


0.0100 


0.100 


(2.60 ±0.76) x 10 17 


chaotic 


0.0100 


0.010 


(1.76 ±0.82) x 10 16 


chaotic+regular 


0.0100 


0.010 


(1.38 ±0.62) x 10 16 


regular 


0.0100 


0.010 


(1.83 ±0.84) x 10 16 


chaotic 


0.0100 


0.001 


(1.07 ±0.57) x 10 16 


chaotic±regular 


0.0100 


0.001 


(8.69 ± 3.80) x 10 15 


regular 


0.0100 


0.001 


(1.56 ±0.65) x 10 16 


chaotic 



TABLE IV. Maximum power emitted as a gravitationala waves with a = 0.01 for three energies 



E = 0.1,0.01 and 0.001. 
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a 


a Q 


T 


A 




AE/ E max 


AE/ E min 


(Ia) 


0.0001 


1 


8.63 x 10~ 3 


4.99 x 10" 


-16 


4.32 x 10~ 4 


1.21 x 10" 5 


(lb) 


0.0001 


23.7 


1 


1.30 x 10" 


-19 


5.0 x 10~ 5 


3.5 x 10~ 7 


(Ic) 


0.0010 


1 


3.23 x 10~ 2 


1.21 x 10- 


-17 


6.94 x 10" 1 


1.87 x 10" 1 


(Id) 


0.0010 


9.80 


1 


2.80 x 10" 


-20 


4.98 x 10~ 2 


1.34 x 10~ 2 


(Ie) 


0.0100 


1 


8.63 x 10~ 2 


4.99 x 10" 


-19 


1.12 x 10" 1 


4.31 x 10" 1 


(If) 


0.0100 


5.12 


1 


8.39 x 10" 


-21 


2.18 x 10~ 2 


8.39 x 10~ 2 



TABLE V. Some set of values for a , T and A when typical size correspond to our Solar System 



with ai = 


1.5xlO n 


m and T 2 


= 3.1 xlO 7 s (see 


text). In the last two columns, we uses E max = 0.10 


and E m i n 


= 0.001 














a 


a Q 


To 


A 


AE/ E max 


AE/ E rnin 


(Ila) 


0.0001 


1 


1.47 x 10~ 2 


2.42 x 10- 11 


3.56 x 10 1 


9.67 x 10- 1 


(lib) 


0.0001 


16.7 


1 


2.12 x 10~ 14 


2.12 


5.77 x 10" 2 


(lie) 


0.0010 


1 


5.51 x 10~ 2 


3.31 x 10~ 13 


3.25 x 10 3 


8.76 x 10 3 


(Hd) 


0.0010 


6.91 


1 


2.63 x 10~ 15 


4.68 x 10 3 


1.26 x 10 3 


(He) 


0.0100 


1 


0.14 


2.45 x 10" 14 


4.75 x 10 5 


1.7 x 10 6 


(Ilf) 


0.0100 


3.59 


1 


9.92 x 10- 16 


2.6 x 10 3 


1.0 x 10 4 



TABLE VI. Some set of values for a a , T Q and A when typical size correspond to a bynary pulsar 
with a2 = 2.0 x 10 9 m and T 2 = 2.8 x 10 4 s (see text). In the last two columns, we uses E max = 0.10 
and E m i n = 0.001 
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a=10" 2 E = 1(T 1 




Figure 1 

"Gravitational Waves from The Newtonian plus Henon-Heiles System" 
by Fernando Kokubun 




Figure 2 

"Gravitational Waves from The Newtonian plus Henon-Heyles System" 
by Fernando Kokubun 



a=10" 2 E = 10" 3 




Figure 3 

"Gravitational Waves from The Newtonian plus Henon-Heiles System" 
by Fernando Kokubun 



